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A detailed analysis of muon-spin rotation (/iSR) spectra in the vortex state of type-II supercon- 
ductors using different theoretical models is presented. Analytical approximations of the London 
and Ginzburg-Landau (GL) models, as well as an exact solution of the GL model were used. The 
limits of the validity of these models and the reliability to extract parameters such as the magnetic 
penetration depth A and the coherence length ^ from the experimental /iSR spectra were inves- 
tigated. The analysis of the simulated /iSR spectra showed that at high magnetic fields there is 
a strong correlation between obtained A and ^ for any value of the Ginzburg-Landau parameter 
K = A/^. The smaller the applied magnetic field is, the smaller is the possibility to find the correct 
value of ^. A simultaneous determination of A and ^ without any restrictions is very problem- 
atic, independent of the model used to describe the vortex state. It was found that for extreme 
type-II superconductors and low magnetic fields, the fitted value of A is practically independent of 5. 
The second-moment method frequently used to analyze fiSH spectra by means of a multi-component 
Gaussian fit, generally yields reliable values of A in the whole range of applied fields Hd <^ H ^ 11^2 
{Hci and 11^2 are the first and second critical fields, respectively). These results are also relevant 
for the interpretation of small-angle neutron scattering (SANS) experiments of the vortex state in 
type-II superconductors. 



I. INTRODUCTION 

The muon-spin rotation (/iSR) technique is one of the 
most powerful and unique tool to study the internal mag- 
netic field distribution P{B) associated with the vortex 
lattice in type-II superconductors (see, e.g. Refs. |T1[21|3|. 
In the vortex state for an applied magnetic field H > Hd, 
or B > {Hci and B are the first critical field and the 
magnetic induction in a sample, respectively) [4J the en- 
ergy of the surface, separating normal and superconduct- 
ing fractions of the sample becomes negative and the field 
penetrates the sample in the form of quantized flux lines, 
called vortices each of them containing an elementary flux 
quantum ($0 = h/2e ~ 2.0678 • IQ-^^ Wb) 0. In the 
case of small pinning these vortices arrange themselves 
in a regular vortex lattice called a flux-line lattice (FLL) 
[5]. The distribution of the internal magnetic fields P{B) 
inside the superconducting sample in the vortex state is 
uniquely determined by two characteristic lengths, the 
magnetic field penetration depth A and the coherence 
length ^. From /iSR experiments P{B) profiles are ob- 
tained by performing a Fourier transformation of the /iSR 
time spectra. There are different approaches to ana- 
lyze fxSR data. Generally, the magnetic field penetration 
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depth A is determined from the second moment (AB^) of 
the internal field distribution P{B) [6l H [HI El UHl [HI IIS] ■ 

For an isotropic extreme type-II superconductor (A ^ f) 
it was shown that (AB^) oc A"'' [H]. The more advanced 
approaches that allow to obtain not only A, but also the 
coherence length ^, require a theoretical model for the 
spatial variation of the internal magnetic field B{r) (r 
is spatial coordinate). An essential requirement of the 
model is that it must account for the finite size of the 
vortex cores. So far, the internal magnetic field distribu- 
tion P{B) measured by /iSR was analyzed assuming an- 
alytical models for B{v) based on London and Ginzburg- 
Landau (GL) theories. The London theory provides the 
simplest approach to model the FLL. Since London the- 
ory does not account for the finite size of the vortex cores, 
a cutoff factor derived from GL theory must be inserted 
into the analytical London expression for B(v) to correct 
for the divergence of i3(r) in the vortex core |3J[T3]. The 
GL theory has the spatial dependence of the order pa- 
rameter build in and thus provides a phenomenological 
description of the magnetic field profile in the vortex core 
region. Abrikosov [5J predicted the vortex core state from 
a periodic solution of the GL equations near the second 
critical field Bc2 = Ii-oHc2 and provided an approximate 
analytical solution of these equations for an isolated vor- 
tex for fields of the order of Hd- Clem [15^ proposed a 
variational method to solve the GL equations that was 
further extended by Hao et al. [16' . A simplified version 
of this model for X/£^^ 1 was developed by Yaouanc et 
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al. [13], and is often used in the literature [5]. 

The London and the GL models were widely applied 
to determine values of A and ^ from measured /iSR time 
spectra taken in the mixed state of type-II superconduc- 
tors [iiiiiniiiaiiiiiniisiiiiiiiaiMiiisiiisiisziiiHi 

[501 [SB • We should emphasize, however, that despite 
of the broad usage, the limits of validity of these mod- 
els and the reliability of the parameters extracted from 
the fits are not much discussed in the literature. The 
main purpose of the present paper is to address these 
basic questions. The paper is divided into two parts. In 
the first part we briefiy describe the models often used 
for the analysis of /iSR spectra: The London model with 
Gaussian cutoff (LG model) , the modified London model 
(ML model) , and the analytical Ginzburg-Landau model 
( AGL model) . These models are compared with the most 
precise model based on the iterative method for solv- 
ing the Ginzburg-Landau equations developed recently 
by Brandt the so called numerical Ginzburg-Landau 
model (NGL model). P{B) profiles for various sets of A, 
^, and magnetic field B were first simulated by means 
of the NGL model and then analyzed within the frame- 
work of the LG, ML, and AGL model. For further dis- 
cussions, it is convenient to define the reduced magnetic 
field b — B/Bc2- It was found that the ML model can be 
used only for low magnetic fields (& ^ 0.1), while both the 
AGL and the LG model yield reliable results almost in 
the whole magnetic field range. However, the values of A 
and ^ obtained by means of the AGL and the LG model 
deviate systematically from the initial parameters used 
for the simulated P{B) profiles for magnetic fields in the 
range 0.01 ^ b < 1. It was also shown that for b < 0.01 
the P{B) profiles do not depend on the coherence length 
^. In the second part of the paper we present a system- 
atic analysis of simulated /iSR time spectra (with typical 
statistics used in real /zSR experiments) by means of the 
LG model. In the whole field range (0 < 6 < 1) and for 
any values of the Ginzburg-Landau parameter k — X/^ 
there is a strong correlation between the values of A and 
f determined from the fit. This implies that an analysis 
of /iSR data using this approach, without taking into ac- 
count these correlations, may lead to substantial errors 
in the determination of the absolute values of A and ^, 
and even may result in unphysical dependencies of A and 
^ on magnetic field and temperature. In addition, the 
second moment method applied to a multiple Gaussian 
fit was tested in order to check how reliably the pene- 
tration depth A can be determined by this method. In 
particular, the influence of the number of Gaussians used 
in the multi-Gaussian fit on the quality of the fit was in- 
vestigated. For typical statistics used in the experiment 
and practically in the whole field range (0 < & ^ 1), the 
second moment method applied to a multi-Gaussian fit 
may provide correct values for A within a few percent. 

The paper is organized as follows: In Sec. II various 
theoretical models used to analyze /iSR data are briefiy 
described. The dependence of the magnetic field distribu- 
tion P{B) on A, ^, 6, and the Gaussian smearing param- 



eter (Tg, as calculated within the LG model, is discussed 
in Sec. III. In Sec. IV we compare the results obtained 
by means of the models described in Sec. II for the case 
of an extreme type-II superconductor (k = A/^ ^ 1). 
Sec. V comprises the studies of the simulated /iSR data. 
The simulated /iSR spectra were analyzed by means of 
the various models described in Sec. II in order to search 
for possible correlations between the parameters, such as 
A, ^, and ag. The conclusions follow in Sec. VI. 



II. MODELS FOR DATA ANALYSIS 

As mentioned in the introduction, the simplest and 
the most widely used approach for analyzing /iSR data is 
based on the relation between the magnetic penetration 
depth A and the second moment (AB^) of the internal 
field distribution P^d{B) of the ideal FLL |4llT3ll32]: 
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Here, C is the proportionality coefficient depending on 
the value of the reduced magnetic field b — {B)/Bc2 [{B) 
is the first moment of PidiB)] and the Ginzburg-Landau 
parameter k j4l[T3l[32]. In order to estimate (Ai?^) one 
often assumes that Pid{B) is a sum of N Gaussian dis- 
tributions (generally, N = 1,2,3) [551 [Mj: 



P^d{B) 



7p 



N 



A,, 



exp 



B-B, 



(2) 

where Ai, Bi, and 17^/7^ are the weight factor, the first 
moment, and the standard deviation of the i—th Gaus- 
sian component, respectively. 7^ = 27rx 135.5342 MHz/T 
is the muon gyromagnetic ratio. The first and second 
moment of Pid{B) are then readily obtained [551 IM] : 
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With modern computers it became possible to develop 
models that allow to calculate Pid{B) for a FLL as a func- 
tion of various parameters, such as magnetic penetration 
depth, coherence length, applied magnetic field, and FLL 
geometry (rectangular or hexagonal) [51 [H [T71 [55] . The 
London models (with different cut off factors) provide 
the simplest and fastest way to calculate Pid{B) for the 
analysis of /iSR data for k 3> 1 [3j. Better approxima- 
tions of Pid{B) for small values of k and fields closer to 
Bc2 can be obtained by the AGL model [lH [1^ . Strictly 
speaking, Ginzburg-Landau theory is only valid in the 
neighborhood of the phase boundary Tc{B) of a type- 
II superconductor. However, it is generally assumed. 
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that Ginzburg-Landau models are also good approxima- 
tions for any field and temperature. The results obtained 
by the NGL model correspond to the minimum of the 
Ginzburg-Landau free energy, whereas other models de- 
scribed in this paper are just approximations to the NGL 
model. Therefore, the NGL model will be used as a refer- 
ence for comparison with the AGL, ML, and LG models. 
A relatively simple method to calculate Pid{B) within the 
framework of the NGL model was proposed by Brandt 
[11136]. 

In the LG, ML, AGL, and NGL approximation the spa- 
tial distribution of the magnetic field in the mixed state 
of a type-II superconductor is described by the Fourier 
expansion: 



100r 
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(i3)^exp(-*Gr) Bg(A,0. 
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(5) 



Here, r is the vector coordinate in a plane perpendicular 
to the applied field. The origin of the coordinate sys- 
tem is in the center of a vortex core (see, e.g. Ref. 
G = 47r/-\/3a(m'\/3/2, n + m/2) are the reciprocal lattice 
vectors for the hexagonal FLL, a is intervortex distance, 
Bq, are the Fourier components, and to, n are integer 
numbers. 

For the LG model the Fourier components Bg are [51l35j: 
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For the ML model the Fourier components Bq, are given 

by mm- 
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For the AGL model the Fourier components Be, are [T¥', 
[16]: 
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S (f2^ + x^G^)y^K,{^-^U' 
where foo — ^ ~ b'^, and 

e„ = aV2 - — )(i + b^y/'[i - 2&(i - bf]^/\ 

K 

Here, Ki{x) is the modified Bessel function. For applied 
magnetic fields H ^ Hd the relation hqH ~ (B) holds 
t36j. Finally, for the NGL model no analytical solution 
for the Fourier components Bg exists. They are deter- 
mined numerically jU [33] . 

From the known spatial distribution of the magnetic 
field B{r) in the mixed state one can extract the internal 
magnetic field distribution Pid{B) for the ideal FLL by 
means of the following equation: 



P^d{B) 



J S{B - B')dA{B') 



(9) 




FIG. 1: Example of a spatial distribution of the magnetic field 
B{r) and the corresponding local magnetic field distribution 
Pid{B) for an ideal hexagonal FLL determined by the NGL 
method. The parameters used for the calculations are: A = 
50 nm, 5 = 20 nm, and (B) = 0.35^2 ^ 246.8 mT, and 
intervortex distance a = 69.5 nm. 



where dA{B') is the elementary area of the FLL with a 
field B' inside, and the integration is over a quarter of 
the FLL unit cell [13] • In order to take into account pos- 
sible random deviations of the flux core positions from 
their ideal ones (vortex disorder) and/or possible broad- 
ening of the /iSR spectra due to nuclear depolarization, 
one may convolute the ideal distribution Pid{B) with a 
Gaussian distribution 1351: 



P{B)^-^ [ PUB')exp 



B- B' 



dB' 



\J^liU g 

(10) 

where Og is the width of the Gaussian distribution. The 
relation between dg, vortex disorder, and nuclear depo- 
larization is described in Sec. IIIID] 

The /iSR time spectra can be further simulated by per- 
forming the Fourier transform of P{B) convoluted with 
the Gaussian function given in Eq. ([To|: 



Ae''*' J P{B)e'^'^^'dB , 



(11) 



where A and are the initial asymmetry and the phase of 
the fiSH time spectra, respectively. For the calculations 
of the spatial magnetic field distribution B{r) in the FLL 
31x31 Fourier components of the magnetic field and the 
reciprocal vector G were used. This allows to calculate 
the second moment of P{B) with a precision of better 
than 10~^. The integral in Eq. ([9| was calculated nu- 
merically over a quarter of the FLL unit cell, divided in 
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approximately 100x100 equal pixels, depending on the 
mean magnetic field {B) (see Fig. 1 of Ref . i29i) . 

Equations (|5] - [s]) are only valid for isotropic super- 
conductors or superconductors with axial symmetry with 
the external magnetic field applied along the symmetry 
axis. In the present study we mostly concentrate on the 
case of extreme type-II superconductors (k ^ 1), such 
as the cuprate high-temperature superconductors. Since 
the qualitative behavior of Pid{B) as a function of var- 
ious parameters is essentially the same for a hexagonal 
and a square FFL, we will consider here only the hexag- 
onal case. 

In Fig. [l] we plot the spatial distribution of the mag- 
netic field B{v) in the mixed state and the corresponding 
local magnetic field distribution Pid{B) for A = 50 nm, 
^ = 20 nm, and {B) = 0.35^2 ^ 246.8 mT, as determined 
by the NGL model. The ideal FLL has three characteris- 
tic fields: (i) the maximal field B^^ax corresponds to the 
field in the vortex core, (ii) the field at the peak of Pid{B) 
is the saddle point field Bsad (located in the middle be- 
tween neighboring vortices), and (iii) the minimal field 
Bmin is in the center of the triangle of vortices forming 
the hexagonal FLL [H] . Instead of the full local magnetic 
field distribution Pid{B) we will use these characteristic 
fields to discuss the dependence of the shape of Pid{B) 
on different parameters. 



III. DEPENDENCE OF P{B) ON A, (B), AND 



In this section we concentrate on the analysis of the 
shape of P{B) given in Eq. ^ as a function of penetra- 
tion depth A (Sec. [Ill A I, coherence length ^ (Sec. IIIBl 



smearing width cjg (Sec. HID I 



mean magnetic field {B) (Sec. IIICl, and Gaussian 



A. Dependence of Pid{B) on A 

In Fig. [2] we show examples of the magnetic field dis- 
tribution Pid{B) for different values of the magnetic pen- 
etration depth A at constant mean field {B) — 0.3i?c2 = 
246.8 mT and coherence length ^ = 20 nm, as calculated 
by the NGL model. The region between the minimal and 
the mean field {B) is most important, because the high 
field tail is usually below the noise level of experimental 
/iSR spectra and is generally not observed, especially at 
low fields and for k ^ 1. Our calculations show that 
the differences between the characteristic fields and the 
mean field {B) are proportional to 1/A^. This is in full 
agreement with the results of Sidorenko et al. |37j who 
obtained for applied fields Hd ^ iJ <C i?c2 and k 1 
(in this case {B) ~ fJ-oH) in the London approximation 
the following expressions: 



SB^ad = B,ad - {B) = - -($o/4^A2) In 2 , 



SB^ax = B,nax ' (B) - 2($o/47rA2) In 



(13) 
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Here, a is the intervortex distance, and K = K{l/^/3) ~ 
1.926 is the complete elliptic integral of the first kind [57] . 
Hereafter, for convenience the quantities SBmim 5Bsad, 
and SBjnax defined above are denoted as characteristic 
fields as well. From Fig. [2] and the pronounced depen- 
dence of the characteristic fields on 1/A^ it is evident that 
the /iSR time spectra strongly depend on A. Therefore, 
it should be possible to extract reliable values of A from 
experimental /iSR data. 
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FIG. 2: Local magnetic field distribution Pid{B) for an ideal 
hexagonal FLL obtained by the NGL model for different val- 
ues of A, at fixed ^ and applied field Bapp ~ (5). The curves 
are normahzed, so that J Pid{B) dB = 1. Note that the shape 
of Pid{B) strongly depends on A. 



B. Dependence of Pid{B) on ^ 

Figure [3] shows the ^ dependence of the characteristic 
fields SBa {a — min, sad, max) normalized to 4>o/47rA^ 
[cf. Eqs. p^pl| ] for a set of different mean fields (B), 
as obtained by the LG model. All the characteristic 
fields SBa disappear at ^ > (<i>o/27r(i3))^/^ ($0 is the 
fiux quantum), where superconductivity vanishes. Below 
a certain value of ^ the characteristic fields SBmm and 
SBsad are independent of whereas 5Bmax still depends 
on ^. However, in real /iSR experiments SBmax cannot 
be determined out of the noise level at low (B) . There- 
fore, at these low values of ^ /iSR spectra are practically 
independent of ^. In order to get a feeling what this 
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FIG. 3: ^ dependence of the characteristic fields SBa {a = 
min, sad, max) normahzed to <I>o/47rA'^ for a set of different 
apphed magnetic fields {Bapp ~ {B) = O.Of T, O.f T, f T, 
fO T) as obtained by the LG model. Note that there is a crit- 
ical value of ^ below which SBmin and SBsad are practically 
independent of 5. This critical value depends on (B). 




b = 27r<B>^7<D 







FIG. 4: Characteristic fields SBa (a = min, sad, max) 
of Fig. |3] plotted as a function of the reduced field b = 
{B)/Bc2(0 iBc2{0 = <&o/27re') at (B) = 0.1, 1, and 10 T 
(black solid line), and at 0.01 T (blue dotted line). Note that 
all the curves SB^iO oi Fig. [3] at (B) = 0.1, 1, 10 T merge 
to single curves SBa{b). 



means for cuprate superconductors we assume ^ ~ 3 nm, 
a typical value of ^ for cuprates below Tc/2. In this 
case the shape of Pid{B) is almost independent of ^ for 
fields {B) < 0.3 T, where 6B„iin and SBsad saturate (see 
Fig. |3|. It is thus problematic to find the correct value 
of ^ at low magnetic fields. At higher fields the shape of 
Pid{B) strongly depends on ^. Note that in Fig. |3] the 
curves corresponding to the smallest field {{B) = 0.01 T) 
exhibit slightly smaller saturation values of the charac- 
teristic fields SBjnin and SBsad than for the higher fields. 
The reason for this is discussed in Sec. IIII CI 

As shown by Brandt |3] the ideal internal field dis- 
tribution Pid{B) may be expressed by normalized param- 
eters, depending only on k = A/^ and b — {B)/Bc2- In 
a similar way we can plot the curves in Fig. [3] not as a 
function of but as a function of 6 = {B)/Bc2, where 
Bc2{0 — *I'o/27rf^ (the relation obtained from Ginzburg- 
Landau theory). This plot is shown in Fig. |4) All the 
curves of Fig. |3] except the one for the smallest field 
{B) = 0.01 T, fall on the same curve. 



C. Field dependence of Pid{B) 

Before we discuss the field dependence of the charac- 
teristic fields SBa {a — min, sad, max) on various pa- 
rameters, it is useful to define the minimal value of the 
reduced field bmin = Bd/ Bc2 — lnK/2K^ which is needed 
to form a regular FLL. This field corresponds to the limit 
below which the vortices can be considered as well sepa- 
rated and noninteracting. 

Figure [5] shows SB^ (normalized to <i>o/47rA^) as func- 
tion of the reduced magnetic field b for different values of 
Av, as calculated by the LG model. The arrows at SBsad 
correspond to bmin{i^)- This figure looks very similar to 
Fig.|4]and represents actually its generalization. It shows 
how SBa depends on all three parameters A, ^, {B), and 
not only SBa as a function of {B). Since Fig. |5] demon- 
strates the dependence of SBa on all the parameters it is 
the basis of all further discussions. At high values of 6, 
the characteristic fields SBa{{B)) and SBa{^^) coincide, 
but at lower fields they deviate substantially (dependence 
of SBa on a parameter x means that all other parame- 
ters, except X, are fixed). The reason for this is obvious. 
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FIG. 5: Characteristic fields 5Ba {a — min, sad, max) cal- 
culated by the LG model as a function of the reduced field b 
(dashed, dotted and dash-dotted lines) for different values of 
K. The black solid lines represent the curves SBa (C) shown in 
Fig. |4] for {B) = 0.1 T, 1 T, and 10 T. The arrows indicate 
the values of bmin at which (B) — Bd- 



For 5 ^ at constant f or Bc2 the intervortex distance a 
increases, and {B), 6Bmin{b), and 6Bsad{b) tend to zero 
as well. This is the reason for the smaller saturation val- 
ues of SBa at the lowest field (B) = 0.01 T in Fig. |3] 
However, in the case of the ^ dependence, when & ^ at 
constant field (B), only the vortex core size is reduced, 
and the intervortex distance a does not change. This has 
not much infiuence on the internal magnetic field distri- 
bution Pid{B) for K 1 (see Fig. [5]). When k is reduced, 
the deviation of 5Ba{{B)) from 6Ba{{B))i^=oo starts at 
higher values of b. For small values of k the characteristic 
fields SBa{{B)) even do not reach saturated values as in 
the case of high k and small b. Despite of the similarity 
of 6Ba{{B)) and SBa{£,^), the mean field (B) can easily 
be extracted from the fit (unlike ^), since it defines the 
oscillation frequency of the fiSR time spectrum. 
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FIG. 6: Change of the local magnetic field distribution 
Pid(B) for an ideal FLL after convolution with a Gaussian dis- 
tribution of various widths (jg . The following parameters were 
used to generate Pid{B) with the LG model: A = 200 nm, 
f = 4 nm, and {B) = 100 mT. 



by several factors, which generally lead to an additional 
broadening of Pid{B) : (i) The FLL is never ideal, but dis- 
ordered by random pinning effects of the vortex cores, (ii) 
For powder samples of anisotropic superconductors - such 
as the layered cuprate superconductors - the grains usu- 
ally have random shapes, and therefore have anisotropic 
superconducting properties and demagnetization effects 
play a role |33j . (iii) The sample may contain magnetic 
nuclear moments or paramagnetic impurities. Vortex 
disorder and nuclear broadening can be taken into ac- 
count by convoluting the ideal internal field distribution 
Pid{B) with a Gaussian distribution of width (see Eq. 10 1 
[ni[351[3H]: 



VD 



(15) 



where avD and ctat are the contributions to the Gaus- 
sian broadening of Pid{B) due to vortex disorder and 
nuclear depolarization, respectively |45j . For k 3> 1 the 
standard deviation of the vortex core positions from the 

1/2 

ideal positions in the FLL (s^) is related to (Jvd by 
the following equation [T7| : 



uvD « A ^ • (1 - 6) 



1/2 



(16) 



D. Dependence of P{B) on ag 

In reality the internal magnetic field distribution in 
the mixed state of a type-II superconductor is infiuenced 



Here b — {B)/Bc2, and a is the intervortex distance. 

Figure [6] shows examples of P{B) for A = 200 nm, 
^ = 4 nm, (B) = 0.1 T, and for various values of ag, 
calculated by means of the LG model. It is obvious that 
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FIG. 7: Characteristic fields SBa (a = min, sad, max) as a 
function of the reduced field b — 27r(i3)^^/"l?o for a = oo, as 
calculated by the LG, ML, AGL, and NGL model. 



with increasing degree of disorder the Van Hove singu- 
larities in the ideal internal field distribution Pid{B) are 
smeared out. Note that the low-field part of Pid{B) is 
mainly truncated by the Gaussian smearing, whereas the 
high-field tail is nearly not affected. 



IV. COMPARISON OF DIFFERENT MODELS 

In this section the different models (LG, ML, AGL, 
and NGL) discussed in this work are compared. For this 
purpose the NGL model, is used as a reference model to 
describe the mixed state of a type-II superconductor. In 
Sec. II we showed that the characteristic fields SBa for 
1 may be represented by single curves (see Fig. |5|. 
Figure [7] shows the characteristic fields 5Ba as a func- 
tion of the reduced magnetic field h = {B)/Bc2 in the 
limit of K ^ oo as calculated by the LG, ML, AGL, and 
NGL model. For small b values 6Bmin{b) and SBsad{b) 
coincide for all models. Deviations of the AGL and LG 
model from the NGL model appear above b « 0.01. Al- 
though the AGL and LG models may well fit the /iSR 
spectra simulated by the NGL model, the fitted values 
of ^ may deviate substantially from the real values for 
reduced magnetic fields b > 10~^. This systematic de- 
viation increases with increasing magnetic field. For the 



LG model, in contrast to the AGL model, the systematic 
errors even change sign with increasing magnetic field. 
One should note, that for k > 5 there is no advantage of 
using the AGL model instead of the LG model. Of all 
the models the ML model approximates best the NGL 
model up to about b ~ 0.1, in agreement with previ- 
ous results of Brandt 13J. However, at higher fields this 
model substantially deviates from the NGL model. The 
ML model was often used to analyze experimental data 
for 6 > 0.1 [21 [29]. We found that fiSK spectra simulated 
by the NGL model in the range = 0.1 — 1 may well be 
fitted by the ML model. But for b > 0.1 the values of ^ 
extracted from the simulated ^SR spectra are artificially 
reduced compared to the real values of ^ (see Fig.[7]). For 
smaller values of k, the characteristic fields SB^ for the 
ML, AGL, and NGL models behave similarly to those 
of the LG model (see Fig. |5|. Namely, for high reduced 
fields b all the curves with different values of k coincide. 
The smaller k is, the higher is the reduced field b when 
they start to deviate from the curves shown in Fig.jTj The 
only exception is the AGL model for which for k < 5 the 
curves become closer to the NGL curves. Here we should 
mention that our results obtained with the NGL model 
are in full agreement with the calculations of Brandt [4j . 



SIMULATION AND FITTING OF /iSR 
SPECTRA 



In order to check the conclusions we made in the pre- 
vious sections "experimentally", /iSR time spectra with 
known parameters ((_B), A, f , Cg) were simulated using 
the LG and the NGL models. For the simulation of 
the ^SR experiment a transverse-field (TF) configura- 
tion with two positron detectors Dl and D2 located 
on opposite sides of the sample was used. The num- 
ber of positrons detected by the detector Dl at time 
ti = i - At is Ni{ti)At {i = 1, 2, 8000; At was chosen 
to be 1.25 ns, corresponding to a typical time-resolution 
for the GPS spectrometer at the Paul Scherrer Institute, 
Switzerland) . This positron count number obeys Poisson 
statistics, and the standard deviation is simply given by 
y^Ni{ti)At. In the ideal case of no noise, the detector 
Dl would detect the signal ni{ti) = noe~*'/'^f' [1 -I- P{ti)], 
where hq is a constant depending on the number of 
muons detected (statistics) and the time interval At, 
— 2.197019(21) /is is the muon lifetime, and P{ti) is 
the noiseless /iSR time signal [see Eq. (11)] to which 
noise has to be added. The signal monitored by de- 
tector Dl can be simulated by the equation A^i(ti) — 
TT-i^ti) + ni{ti)gi, where gi is a random number gener- 
ator obeying Gaussian statistics with standard deviation 
and variance equal to 1. A similar signal but with op- 
posite phase is registered by detector D2. In analogy to 
real /iSR experiments one can calculate the asymmetry 
[HIS] A{t,) = [Ni{t,)-N2{t,)]/[Ni{t,) + N2{t,)], yielding 
P{ti) with "experimental noise". The /zSR time spec- 
tra were simulated according to the procedure described 
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above with total statistics of 20 million events, a value 
typically used in real experiments. 

The simulated /iSR time spectra were then analyzed 
as follows: 

(1) The /iSR spectra simulated by the NGL model were 
analyzed by the second-moment (SM) method. 

(2) The ^SK spectra simulated by the LG model were 
analyzed by a fitting procedure based on the LG model. 

According to the discussions made in the previous sec- 
tions the following important points emerge: 

(1) How reliable is the second moment obtained by a 
multi-Gaussian fit of the /iSR spectra (see Eq. (17l) and 
the value of the penetration depth A extracted from the 
second moment? 

(2) Is there a correlation between ag and 1 /A^, since both 
of them influence the second moment of the fiSH spec- 
trum? 

(3) Is it possible to extract reliable values of ^ from /iSR 
spectra at low magnetic fields b ~ 10""^? 

(4) Is there a correlation between A and ^ at high fields? 
Since for b « 0.1 — 0.9 both parameters strongly influence 
the characteristic fields SBa- 



A. Test of the second-moment method 

In this section the second-moment (SM) method is 
tested by analyzing ^SR time spectra simulated with the 
NGL model with well defined parameters {{B), A, ^, (Tg). 
The SM method is theoretically well described in the lit- 
erature lU 1131 [37] and was extensively used to extract 
the magnetic penetration depth of extreme type-II su- 
perconductors from fiSR spectra [6l [T) [8] [9l [101 [TT| [T2] . 
In the framework of this method the /xSR time spectra 
are usually fitted to a sum of N Gaussian components 



N 

P{t) = ^ A, exp(-(T^V2) cos{jf,B,t + (j)). (17) 

i=l 

Here is initial phase of the muon beam, and Ai, ai, 
and Bi are the asymmetry, the relaxation rate, and the 
first moment of the z— th Gaussian component, respec- 
tively. From the fit parameters Ai, ai, and Bi one read- 
ily obtains the first and the second moment of P{B) 
from Eqs. ^ and Q, respectively. Using Eq. ([T]) one 
finds the penetration depth A. Here the serious ques- 
tion arises: How reliable is the value of A obtained by 
the SM method using a multi-Gaussian fit? In order to 
investigate this question /iSR time spectra were simu- 
lated by the NGL model for an extreme type-II super- 
conductor (such as the cuprate superconductors) with a 
Ginzburg-Landau parameter k = 50 3> 1. The temper- 
ature dependence of A was assumed to follow the rela- 
tion (two-fluid model): A(T)-VA(0)-2 = [1 - (T/Tc)^] 
with Tc — 22.5 K and A(0) — 200 nm (zero-temperature 
penetration depth). This approximately corresponds to 
the temperature dependence of A in the strong-coupling 



BCS case (3S]. In the first step we assume that there is 
no vortex disorder {avD = 0) and no nuclear depolar- 
ization ((Tjv — 0) present [cf. Eq. (15l]. The simulations 
were performed for three different magnetic fields: (B) ~ 
0.05, 0.5, and 5 T [46 . This corresponds to the reduced 
fields b 0.0024, 0.024, and 0.24, an extremely small, 
an intermediate, and a high magnetic field, respectively 
(see Fig.™. Since Bc2 is decreasing with increasing tem- 
perature (1/^^ oc Bc2', 1/C^ 1/A^ at constant k), the 
analysis of the /j,SR spectra for B — 5T makes sense only 
up to 21 K where i3c2(21i^) ~ (B) and superconductiv- 
ity disappears. Noisy /iSR time spectra were simulated 
with the parameters A, ^, and {B) as described above. 
For the technical parameters the following typical values 
were used: statistics 20 • 10®, asymmetry A = 0.2, and 
phase = 0. 

For the analysis of the simulated fiSH spectra, P(t) in 
Eq. (17 1 was approximated by a sum of N — 1,2,3,4,5 



Gaussians in order to check the reliability of the result 
obtained by a multi-Gaussian fit. The number of Gaus- 
sians N determines the quality of the fit. N should be 
increased from 1 until (normalized to the degrees of 
freedom) is close to 1 within statistical scattering. The 
total asymmetry A — J2i=i ^^"^ the phase (p of P{t) 
in Eq. ( 17 1 were assumed to be known and were fixed 
in the fitting procedure. According to our experience, 
in order to reduce the scattering of the fitted values of 
the second moment, one should fix the asymmetries Ai of 
the individual Gaussians to their average values obtained 
by a fit with all parameters free. From the first and the 
second moments of the individual Gaussians one can cal- 
culate the second moment (Ai?^) of the /iSR spectrum 
using Eq. ^ , that corresponds to the second moment of 
Pid{B) of an ideal FLL. The magnetic penetration depth 
A is readily obtained from (AB^) with Eq. ([l]). The result 
of the analysis of the simulated /iSR time spectra for an 
ideal FLL is shown in Fig. [8| For the 0.05 T data the fi- 
nite value of K = 50 was taken in account in the coefficient 
C in Eq. ^ 0]. One can see that the smaller the field is 
the more Gaussians N are needed to describe the spec- 
tra. Whereas for 5 T = 2, 3 Gaussians are sufficient 
to reproduce the spectra, A^ = 4 and N = 5 are required 
for 0.5 and 0.05 T, respectively. Note that the scatter- 
ing of the data points increases with increasing number 
of Gaussians A^. Although at 0.05 T the fitted values of 
1 /A^ deviate systematically from the real values by a few 
percent, the qualitative behavior of 1/A^(T) is the same. 
As will be shown below by adding a Gaussian smearing 
CTg to the /iSR spectra, the scattering is reduced, and a 
smaller number of Gaussians A^ are needed to describe 
the spectra. Fig. [9] demonstrates how the local magnetic 
field distribution Pid{B) for an ideal FLL can be approx- 
imated by A = 5 Gaussians. Although not all the details 
of Pid{B) are reproduced, the overall agreement is good, 
in particular the second moment. 

In order to test the second moment method under more 
realistic conditions one should add a Gaussian smear- 
ing CTg to the fiSK spectra [cf. Eq. (lOl]. According to 
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FIG. 8: Fit results for obtained by the second moment 
method. The noisy spectra for the three different fields of 
0.05, 0.5, and 5 T were simulated by the NGL method for an 
ideal FLL as described in the text, and then analyzed by a 
multi-Gaussian function with different number of Gaussians 



{N = 1, 2, 3, 4, 5) as defined in Eq. mj. The black solid hues 
correspond to the real values of A""^ used for the simulation. 



Eq. (15 1, we assume for the further discussions that ag 



is composed of two components: ag = (cy^ + cr^)-^/^, 
where (Tvd denotes the temperature dependent smear- 
ing due to vortex disorder, and ctat is the temperature 
independent smearing due to nuclear depolarization [cf. 



1 /2 

Eq. (16)]. For a constant vortex disorder (s^) /a = 



161 



consf7(rigid vortex lattice) and reduced field b, the re 
lation avD oc 1/A^ holds. As is obvious from Eq. ( 
for (B) = 5T the term (1 - 6) = 0.76 substantially c 
ates from unity and has to be taken into account in the 
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FIG. 9: Comparison of the ideal field distribution Pid{B) 
simulated with the NGL model at 0.5 T and 1 K (blue line) 
with P{B) obtained by a Gaussian fit with N — 5 (red dashed 
line). The five individual Gaussian components used for the 
fit are also shown (black dotted lines). The inset shows the 
same plot but on a semi-logarithmic scale. 



simulation of /iSR data. 

For the simulations of the smeared /xSR spectra the fol- 
lowing values for single crystal Lai.83Sro.i7Cu04_5 were 
used: avD = /3/A^ (3 = 2.585 • 10* (mT-nm^) and 
<^N — <^Cu — 0.27 mT HO]. Noisy /^SR time spectra were 
simulated with the parameters A, ^, (i?), and Ug as de- 
scribed above. For the technical parameters the following 
typical values were used: statistics 20 • 10^, asymmetry 
A = 0.2, and phase = 0. 

The total second moment (Ai?^)f of a /iSR spectrum 
with Gaussian smearing Og obtained by multi-Gaussian 
fit contains three components: 



(Ai?2)t 



•VD 



(18) 



where (AB^), Uy^j, and cr^ are the second moments 
due to the internal field variation in the ideal FLL, the 
vortex disorder, and the nuclear depolarization, respec- 
tively. In order to obtain A from the total second mo- 
ment measured in real experiments one determines ct/v 

AB2) in Eq. ([T8|, 

41 



above Tc, and assumes that cTy^, ^ 



(AB2 



(AB2 



•n 



From the first 



and 



the second moments of the individual Gaussians one can 
calculate the total second moment (Ai?^)t of the /xSR 
spectrum using Eq. (Q [Note that in Eq. (Q (Ai?^) has 
to be replaced by (AB^)j for the case Og ^0]. By means 
of Eqs. (18) and ([T]) the magnetic penetration depth A 
then is readily obtained. Figure 10 shows the results for 
the penetration depth obtained by the second moment 
method with = 1, 2, 3, 4 Gaussians. Note that a 
single Gaussian does not give reliable results in agree- 
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FIG. 10: Fit results for obtained by the second moment 
method. The noisy spectra for the three different fields of 
0.05, 0.5, and 5 T were simulated by the NGL method includ- 
ing Gaussian smearing Og as described in the text, and then 
analyzed by a mult i- Gaussian function with different number 
ofGaussians (A'^ = 1, 2, 3, 4), as defined in Eq. (17 1. The black 
solid lines correspond to the real values of A~ used for the 
simulation. 



ment with earlier findings [33! • However, with increas- 
ing number of Gaussians N the quality of the fits sub- 
stantially improves. In order to fit the simulated data 
at 0.05 T at least 3 or better 4 Gaussians are required. 
For = 3 there is a systematic deviation of about 10% 
of A^^ from the real value (or 5% for A), whereas for 
= 4 the values of A~^ are scattered within a few per- 
cent around the real ones. For the data simulated at 
0.5 T even = 3 Gaussians are sufficient to describe 
the local magnetic field distribution P{B), and the val- 
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FIG. 11: Comparison of the real field distribution P{B) 
(empty circles) simulated with the NGL model for param- 
eters A = 200 nm, C = 4 nm, {B) = 0.5 T, ag = 0.7 mT (data 
point at r = 1 K for iV = 3 in Fig. [T^) with P{B) obtained 
by a Gaussian fit with N = 3 (red line). The dotted, dashed, 
and dash-dotted lines represent individual Gaussian compo- 
nents used for the fit. It is obvious that a multi-Gaussian fit 
may well describe the real P{B). 



ues of A ^ are systematically shifted only within a few 
percent. Fig. [Tljshows an example of a real internal field 
distribution P{B) {{B) = 0.5 T, T = 1 K of Fig.flOl and 



the reconstructed P{B) obtained from the analysis of the 
simulated ^SR spectrum (2x10'' statistics) using a Gaus- 
sian fit with = 3. It is obvious that three Gaussians 
describe well the shape of the real P{B). The largest 
systematic error in A^^ obtained by a multi-Gaussian 
fit is observed at 5 T (6 = 0.24). At such a high field 
{{B) ~ -Bc2/4) the variation of the internal field is rela- 
tively small (see Fig. [7| , and the Gaussian smearing try £, 
[cf. Eq. (18)] due to vortex disorder becomes essential. 
The second moment of this Gaussian smearing cannot 
be neglected and considerably contributes to the total 
second moment {AB'^)t of P{B). This leads to system- 
atically higher values of A"^ obtained by multi-Gaussian 
fits at high magnetic fields. Note, however, there are 
two reasons why the contribution to the second moment 
due to vortex disorder is reduced with increasing mag- 
netic field, and consequently the systematic error in A~^: 
(1) At high fields avD oc (1-6) [cf. Eq. which 
was taken into account in the simulations of the fiSH 

1/2 

spectra. (2) Vortex disorder, (s^) /a, is expected to 
decrease with increasing magnetic field, because of the 
strong repulsive interaction between the vortices at high 
fields (see e.g. Ref. H]). 

By means of a multi-Gaussian fit it is not possible to 
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FIG. 12: Temperature dependence of 1/A^ determined from 
{AB^), assuming that avo is known as described in the text. 
Triangles: iV = 4, 0.05 T; stars: iV = 4, 0.05 T, corrected for 
finite k = 50 (see text for explanation); circles: A'^ = 3, 0.5 T; 
squares: A'' = 3, 5 T. 



separate (Ai?^) and avo from the measured total second 
moment {AB^)t [cf. Eq. (18l]. Assuming that avD = 
yields a lower limit for A (upper limit for 1/A^), as clearly 
demonstrated in Figs. 10 'b) and 10 ^c) where the values 
of 1/A^ are systematically too large. It is interesting to 
investigate what are the values of 1/A^ after correction 
with the real value of ayo- For this purpose we write 
Eq. (16) in the form avD — /5A^^ and with the help of 
Eqs. ipl and ( 18 ) we obtain: 



(19) 



where e = Cf3^ describes the correction d ue t o vortex 
disorder. The values of 1/A^ plotted in Fig. 10 were ob- 
tained with e = (no vortex disorder correction) . Fig. 12 



shows some of the results of Fig. [TO] after correcting the 
values of 1/A^ with the values of €t=o = 0.050, 0.073, 
and 0.155 for 0.05 T, 0.5 T, and 5 T, respectively (e is 
temperature dependent, since C(6(T),k) is temperature 
dependent). The corrected values of 1/A^ are in good 



agreement with the real values (solid line in Fig. 12 1, ex- 



cept for the data at 0.05 T where a systematic deviation 
of about 5-10% is observed. One of the reasons for this 
deviation is that k = 50 used in the simulations is not 
infinite. This implies that at 0.05 T (6 ~ 0.0024) the pa- 
rameter C~^/^ in Eqs. \l\ and (19) is about 5% smaller 



[4]. The stars in Fig. |12|represent the corrected values of 
1/A^ at 0.05 T, which are only about 3% systematically 
lower than the real values. 



B. Test of the London model with Gaussian cut off 

(LG) 

In order to test the reliability of the advanced methods 
described in Sec. II we first simulated noisy /zSR spectra 
and then fitted them in a similar way as for the sec- 
ond moment method. To avoid systematic errors in the 
fit results it is important to analyze the data with the 
same model as they were simulated. Here we present 
results to the LG model, since it can well approximate 
experimental data in the whole field range (see Fig. [7| . 
Similar results are also obtained with all other models. 
The temperature dependence of the penetration depth 
A was assumed to follow the relation A~^(r)/A~^(0) = 
[1 - (T/Tc)2] with A(0) = 200 nm and = 22.5 K. 
The Ginzburg-Landau parameter n — 50, and Gaussian 
smearing ag{T) = (^^.^(r) + a%y/^ [avo « 1/A2(T)] 
was chosen to be the same as in Sec. |V A[ Again the 
/iSR spectra were simulated for three different mean fields 
(B) = 0.05, 0.5, and 5 T (& = 0.0024, 0.024, 0.24). As 
before the mean field was assumed to be temperature in- 
dependent. Statistics, asymmetry, and phase of the /iSR 
time spectra were chosen to be 20 • 10^, 0.2, and 0, re- 
spectively. 

The results of the fits of the simulated fiSK spectra 
with all the parameters free (except A and (j>) are shown 
in Fig. |13[ Phase (j> and asymmetry A were assumed to 
be known and were fixed to their real values. In order to 
exclude any artificial infiuence on the fitting procedure 
we performed the fits in automatic mode. This means 
that with increasing temperature the fit results for tem- 
perature Ti were used as initial parameters for the next 
temperature T^+i. For Ti (lowest temperature) the cor- 
rect initial parameters were used. As shown in Fig. |13| for 
A~^ we got good results at low fields, however there are 
substantial systematic deviations for k. At the highest 
field there are substantial systematic deviations of the 
fitted values for both A and k, although the goodness of 
fit for these fits is comparable to those for the correct 
parameters (Fig. 13). The values of weighted and 
normalized to the degrees of freedom (w 8000) scatter 
around 1 within 0.04, as expected f or th e present degrees 
of freedom. The black stars in Fig. 13 denote for the 
true values of the parameters for comparison. From all 
the fitted parameters, the values obtained for k deviate 
mostly from the real ones for all the fields. In order to 
check the reason for this strong deviation for k at 0.05 T, 
the fit was performed using different initial parameters. 
For A, (B), and Cg correct values within a few percent 
were obtained, whereas ^ was found to be in the range 
2-13 nm with a very good value of ■ The dependence of 
on ^ for different fields and statistics is demonstrated 
in Fig. [14] It is evident that at low fields and 20 million 
statistics the quality of the fit is practically independent 
of ^ in a very broad range. For the 5 T data not only 
K but also A substantially deviates from the real value 
(see Fig. [T3| . The good agreement at low temperatures 
is misleading, since it is only due to the correct initial pa- 
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FIG. 13: Summary of the fit results for the /iSR spec- 
tra simulated by the LG model. The fitted values 1/A^, k, 
A{B)/{B), and ag, as well as s^re plotted as a function 
of temperature T for three different fields 0.05 T (circles), 
0.5 T (up-triangles), and 5 T (down-triangles) for 20 mil- 
lion statistics. The parameters 1/A^, k, A{B)/{B), and ag 
were free during the fitting procedure. The solid lines cor- 
respond to the true values of the parameters. A{B)/{B) = 
{{B) fit — (B) real) /{B) real dcuotes the relative deviation of 
the fitted value (B) fu from the real {B)reai- For comparison 
was also determined for the real values of the parameters 
(black stars). The black empty diamonds in (a) show a pos- 
sible fit result for at 5 T with extremely wrong initial 
parameters. Note that the error bars of l/A'^, as calculated 
by the fitting program (the function "fit" of the MATLAB 
program was used), are within the point size (~ 2%). 



rameters we set for the lowest temperature. The empty 
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FIG. 14: a function of ^ for the /iSR spectra simulated 

by the LG model with 20 million counts, (B) = 0.05, 0.2, 
0.5, and 5 T, A = 200 nm, ^ = 4 nm, and ag = 0.7 mT. 
For comparison X^iO also shown for the higher statistics 
of 2000 millions at (B) = 0.05 and 0.2 T. The statistically 
scattered minimal value of x^ was normalized to 1. Note that 
the dependence of x^ on ^ at low fields is weak. 



black diamonds in Fig. 13 a) shows the fit result for in- 
tentionally extremely wrong starting parameters. In real 
measurements one never knows the optimal starting pa- 
rameters. We performed fits of the 5 T data at T = 1 K 
with different starting values of A, ^, {B), and ag. It was 
found that for (B) and ag one obtains values close to 
the real ones, however, for A and ^ this is not the case. 
Figure [15] shows the variation of the values of A and ^ 
during the fitting process. The starting values of A and 
^ are indicated by number The numbers [2], [3~| , 
indicate the values of A and ^ after each 5 fit- 
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ting iterations. The maximal number of fitting iterations 
was not restricted. However, 50 iterations were usually 
sufficient, and the fit was terminated when the relative 
changes of all the parameters during the iteration were 
less than 10^^. The final results of the fit eventually 
correspond to local minima of x^. The fit converges on 
a certain A — A(f) curve denoted by the black line in 
Fig. |15| indicating a possible correlation between A and 
^ at high fields. Therefore, we can conclude, that the 
determination of a reliable value of ^ is problematic at 
low fields as expected (see Fig. [S]). At higher fields not 
only the value of ^, but also the value of A may systemat- 
ically deviate from the real value. However, at low fields 
reliable values of (B) , ag , and A may presumably be well 
determined from fits with the advanced models. 

The next step for improving the fitting procedure is 
to restrict some parameters. Based on the results ob- 
tained from the free parameters fits we conclude that 
a good candidate for restriction is ^, especially at low 
fields. One can fix k in order to relate ^ to A via ^ = A/k 
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FIG. 15: Visualization of the fitting process of tlie simulated 
fiSK spectrum for the parameters A = 200 nm, ^ = 4 nm, 
(B) = 5 T, ag = 0.7 mT using the LG model. Number [l] 
indicates the starting values for A and The numbers [2^ , 
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denote the values of A and ^ after each 5 iterations 
of the fitting process. The fit was terminated when relative 
changes of all the parameters were less than 10"''. The solid 
line shows the points in the A vs. ^-plane where the fit is 
finally converging. Note that the fit results for the other pa- 
rameters (B) and cTg are close to the real ones, independent 
of the starting values. 



and let the other parameters A, (B), and ag free. This 
is also reasonable from a theoretical point of view. In 
the BSC approximation k changes not substantially with 
temperature in the weak coupling limit [42^ . The results 
of the fits with the only restricted parameter k (i.e., ^ 
A/k) are shown in Fig. 



16 It is 



was calculated with ^ 
obvious, that the fits are excellent and all the parameters 
are very close to the real ones with only small statistical 
scattering. 



Correlation between a„ and A for small values 



of b 



As shown in Fig. [6] with increasing Gaussian smear- 
ing (Tg the characteristic fields of the internal magnetic 
field distribution Pi^^B) of the ideal FLL are gradually 
washed out, and P{B) tends to become an asymmetric 
Gaussian-like distribution. Therefore, one expects some 
correlation between ag and the inverse square of the pen- 
etration depth 1/A^, since both of them influence the 
second moment of the ^SR spectrum. In order to show 
the possibility of extracting the real values of A and (Tg 
from /iSR spectra, we simulated /iSR data by the ML 
model and then calculated the goodness of fit as a 
function of A and ag with the other parameters fixed 
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FIG. 16: Summary of the fit results for the /iSR spectra sim- 
ulated by the LG model. The fitted values 1/A^ A{B)/{B), 
and (Tg, as well as a^re plotted as a function of tempera- 
ture T for three diflferent fields 0.05 T (circles), 0.5 T (up- 
triangles), and 5 T (down-triangles) for 20 million statistics. 
The parameters 1/A^, A{B) / {B) , and ag were free during the 
fitting procedure, whereas n, was fixed at the real value. The 
solid lines correspond to the true values of the parameters. 
A{B)/{B) = ((B) fu - (B) real) /{B) real describes the relative 
deviation of the fitted value {B)fit from the real one {B)reai- 
For comparison was also determined for the real values of 
the parameters (black stars). 



to their true values. Figure 17 ^a) shows a func- 

tion of A with the other parameters fixed to their true 
values. For the simulated data the following parame- 
ters were used: Aq = 200 nm, = 1 nm, bo = 10^^, 
ago = 0.05, 0.75, 1.5 mT, and 20 million counts. It is 
evident that with increasing ago the error of A extracted 
from the fit increases. Fig. 17 'b) shows a contour plot 
of as a function of A and ag calculated for data with 
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FIG. 17: (a) a-s a function of A with the other parameters 
set to the true values. The following parameters were used for 
the ^SR data simulation: Ao ~ 200 nm, = 1 nm, bo = 10"'', 
20 millions statistics, and a^o = 0.05, 0.75, and 1.5 mT. (b) 
Contour plot of ^ ^ function of A and ag for the data 
simulated with Ao = 200 nm, = 1 nm, bo = 10~^, ago = 
0.75, and 20 million counts (red point). In both figures the 
ML model was used for the calculations, and the statistically 
scattered minimal value of the x^ ^a-s normalized to 1. The 
real values of A and ^ are indicated by the red point. 



Ao = 200 nm, = 1 nm, bo = 10"'^, ago = 0.75, and 
20 millions statistics. This approximately corresponds to 
the case we analyzed before. From the figure we conclude 
that and ag are slightly correlated, but it is possible 
to extract both of them simultaneously if ^ is fixed. This 
agrees well with the results of the analysis performed in 
Sec. EH 



D. Correlation between ^ and A 

For low magnetic fields the dependence of the /iSR 
spectrum on the coherence length ^ is very weak (see 
Fig. 14 1. But with increasing field towards Bc2 the shape 



of the spectrum becomes dependent not only on the pen- 
etration depth A, but also on ^ as well (see Fig. [t]). An 
increase of A and/or ^ causes a decrease of the second 
moment and the characteristic fields. Therefore, it is ex- 
pected that a decrease of A is correlated with an increase 
of ^ in the fitting procedure and vice versa. So far, to our 
knowledge this problem was discussed previously only by 
Riseman et al. |17j . Here we study this problem in more 
detail. For this procedure we determined for simulated 
^SR data as a function A and <^ at fixed {B)q = 0.5Bc2 
and Gaussian smearing = 0.5 mT. We have chosen 
the case of a relatively small kq = 2.5 (Aq — 50 nm and 
^0 = 20 nm) for the data simulation, since the relative 
volume of the vortex cores is large, and therefore it is 
easier to extract ^ from the fits. As before the statistics 
was 20 • 10^. The result of the analysis with the NGL 
model is shown in Fig. 18 where a contour plot of 



as a function of A and ^ with the other parameters fixed 
is displayed. There is indeed a strong correlation be- 
tween A and For A(^) « 58. 68-^2. 14^-0. 127^^ where 
~ 1 is minimal, the fits converge after a few hundreds 
of iterations. Tests showed that for different starting pa- 
rameters the fits were converging in the correlated region 
of ^ = 14 — 24 nm and A « 36 — 63 nm. This region lies 
within the region of x^ < 1-05 (see Fig.n8|. It should be 



noted that for such a high reduced field = 0.5 as used 
here for the analysis, the qualitative dependence of the 
characteristic fields 6Ba on A, ^, and {B) are indepen- 
dent of K = A/^ (see Fig. |5|. Therefore the qualitative 
behaviour of the contour plot of {'^ / -^o , ^ / £,o) in Fig. 18 
is independent of kq = Aq/^o for any kq > 5. 

The situation could be improved for fields b ~ 10~^ 
where the minimal and saddle point fields slightly de- 
pend on ^ (see Fig. [t]), and on the other hand the max- 
imal field is not very large and still depends on ^. How- 
ever, with decreasing field, the vortex core volume sub- 
stantially reduces which is disadvantageous for the data 
analysis as discussed above. Figure 19 shows a simi- 
lar contour plot of X^i-^iO fo'^ noisy data simulated by 
the NGL and ML model with parameters Aq = 200 nm, 
^0 = 4 nm (kq — 50), dgo = 0.7 mT, 20 million statistics, 
and bo = 0.004, 0.02, and 0.1. This corresponds approxi- 
mately to the analysis of the simulated data we discussed 
before. There is a substantial correlation between A and ^ 
at all fields. For 6o = 0.1 the analysis yields ^ = 4(1) nm 
(correlated with A), for bo ~ 0.02 only an upper limit of 
^ ~ 5 nm can be given, and for bo = 0.004 (0.082 T) the 
fit is practically independent of f at 20 million statistics. 
However, at unrealistic high statistics the dependence of 



X on A and f becomes stronger (see Fig. 14 1, and the pre 



cision of the parameters extracted from the fit increases 
as the square root of the statistics. Another way of solv- 
ing this problem was proposed by Riseman ei al. |17j , by 
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FIG. 18: Contour plot of as a function of A and ^ for the 
data simulated with parameters Ao = 50 nm, = 20 nm, 
&o = 0.5, (Jgo = 0.5 mT, and 20 million counts, calculated 
by the NGL model. The fitted parameters A and ^ exhibit 
a strong correlation, which is the reason for the pronounced 
systematic deviations of the fitted values of A~'^ and k at 5 T 
from the real values displayed in Fig. |13| The statistically 
scattered minimal value of is normalized to 1. The real 
values of A and ^ are indicated by the red point. 



simultaneously fitting several spectra measured at differ- 
ent fields with common values of the parameters A and 
^. This has two advantages: (1) It effectively increases 
statistics. (2) Since the correlation curve A(^) depends 
on field (see Fig. 19 1, the total contour graph of X^(''^jC) 
will shrink, allowing a determination of the correct pa- 
rameters. For example in the case of a high value of k, 



and extremely small field b, one can determine the cor- 
rect value of A (independent of the value of and with 
the known value of A it is possible to evaluate a reliable 
value of ^ at high field by means of the correlation curve 
^(A). This procedure can be justified at least for con- 
ventional superconductors. Recently, Landau and Keller 
[32j reanalyzed /xSR data for various conventional super- 
conductors and convincingly demonstrated, that in many 
cases type-II superconductors can be described by a field 
independent penetration depth. The present results are 
also relevant for the interpretation of small-angle neu- 
tron scattering (SANS) experiments in the mixed phase 
of type-II superconductors, since the strong correlation 
between A and ^ is also present in Fourier components of 
the FLL pi. 

We can conclude that in general a simultaneous deter- 
mination of ^ and A from /.tSR spectra without additional 
restrictions is not easy, independent of the model used to 
describe the vortex state. At high k and low fields there 
is practically no dependence of the spectra on ^, and at 
high fields ^ is strongly correlated with A. As demon- 




FIG. 19: Contour plots of x as a function of A and ^ for 
/iSR spectra simulated with the parameters Ao = 200 nm, 
^0 = 4 nm, (Tgo = 0.7 mT, 20 million counts, and 6o = 0.004, 
0.02, and 0.1. The results in (a) were obtained by the NGL 
model, and those in (b) and (c) by the ML model. There is 
a strong correlation of the fitted values of A and ^ for high 
(6o = 0.1) and intermediate (6o = 0.02) fields, but nearly 
no correlation at low fields (6o = 0.004). The statistically 
scattered minimal value of is normalized to 1. The real 
values of A and ^ are indicated by the red point. 



strated in Figs. [3] [S] and [7] this is independent of the 
value of K and the model used. It is important to add 
that in our analysis of the /xSR data by the advanced 
models, we used the same model for the simulations and 
the analysis, what reduces systematic errors to a mini- 
mum. In practice, there is often no adequate model for 
the description of the experimental /iSR spectra as for 
instance for unconventional superconductors, such as the 
cuprate superconductors. In this case additional difficul- 
ties in the data analysis are expected. In the analysis of 
the /iSR spectra in section 5 we did not consider back- 
ground signals arising from impurity fractions/phases in 
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the sample and/or from muons stopping in the sample 
holder or other parts of the spectrometer. These back- 
ground signals may be a hidden source of uncertainties in 
the determination of reliable parameters from fiSH spec- 
tra. The introduction of additional fit parameters in the 
advanced models should be done only with great care, 
since already the existing minimal set of parameters of 
the models are in general difficult to extract. 

VI. CONCLUSIONS 

We performed an analysis of the line shape of /iSR 
spectra of type-II superconductors in the mixed state 
simulated by four different models frequently used: (1) 
the modified London model (ML), (2) the London model 
with Gaussian cut-off (LG), (3) the analytical Ginzburg- 
Landau (AGL), and (4) the numerical Ginzburg-Landau 
(NGL) model. The dependence of the line shape on the 
penetration depth A, the coherence length the applied 
magnetic field B, and the Gaussian smearing parame- 
ter (Tg is in agreement with previous studies jH [35] |37] . 
It is discussed under what conditions these models can 
be used to describe the vortex state in extreme type-II 
superconductors. As a result, the ML model can be ap- 
plied for fields b = B/Bc2 < 0.1 (i?c2 is second critical 
field). On the other hand, the AGL and the LG model 
can be applied in the whole range of fields, but in the 
range b ~ 10~^ — 1 they systematically deviate from the 
NGL model. It was shown that at low fields b < 10^'^ 
there is practically no dependence of the line shape on ^. 
However, with increasing field, there is a strong depen- 
dence of the line shape on both A and but the strong 



correlation between them makes it almost impossible to 
determine A and ^ simultaneously. This is independent of 
K — X/S, and the model used. Additional restrictions for 
^ (or A) are needed to get rid of this correlation for rea- 
sonable statistics. Furthermore, it was shown that it is 
possible to determine A and ag simultaneously, provided 
that ^ is fixed and the correlation between them is not 
too strong. In addition, it was demonstrated that the sec- 
ond moment method (SM), frequently used for /zSR data 
analysis, may yield reliable values for A (within a few per- 
cent) in whole field range < 6 ^ 1, provided that vortex 
lattice disorder is not substantial. A multiple Gaussian 
fit may give reliable values for the second moment and 
may well approximate the local magnetic field distribu- 
tion in a type-II superconductor. In order to substantiate 
these conclusions made above, we performed virtual ex- 
periments by generating noisy fiSH spectra with known 
parameters. The results of a comprehensive analysis of 
these /iSR spectra are in full agreement with the conclu- 
sions made above. 
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